z = 3 antiferromagnetic quantum criticality and emergence of fermionized skyrmions 
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Hertz-Moriya-Millis theory with dynamical critical exponent z = 2 has been proposed to describe 
antiferromagnetic quantum criticality in itinerant electron systems. In this study we show that the 
dynamical critical exponent changes from z — 2 to z = 3 at low temperatures, which results from 
effective long-range interactions between spin fluctuations, generated by Fermi-surface fluctuations 
beyond the Eliashberg framework. We claim that the underlying physics for the z = 3 antiferro- 
magnetic quantum criticality is the emergence of fermionized skyrmion excitations at low energies, 
which form a Fermi surface, referred as a skyrmion liquid, where the interplay between itinerant 
6JQf electrons and skyrmions is argued to allow fermionized skyrmions. We construct a dual field the- 

ory in terms of skyrmion excitations, and obtain the 2 = 3 antiferromagnetic quantum criticality. 
This demonstration suggests the expected but nontrivial consistency between weak-coupling and 
strong-coupling theories. 

oo; 

CN ', PACS numbers: 71.10.Hf, 71.10.-w, 71.27. +a 

! I. INTRODUCTION 

> ■ 

■ At present, it is almost impossible to understand non-perturbative physics based on the weak coupling approach, 
where physically relevant quantum corrections should be taken into account up to an infinite order. Actually, we do 
not know how to describe various strong coupling phenomena such as the Kondo effect, Mott transition, Anderson 
localization, and etc., starting from the Fermi liquid theory and performing the diagrammatic analysis. Dualities 
between weak-coupling and strong-coupling theories have been playing key roles in figuring out non-perturbative 
phenomena, where topologically nontrivial excitations instead of original variables are introduced to describe strong 
coupling physics. Domain wall excitations appear to describe quantum criticality in the dual field theory of Ising 
O ■ model, identified with fermions, where the transformation from the Ising variable to the fermion field is highly 
nonlocal and nonlinear Unfortunately, duality transformations are not well defined in most cases. It is also almost 
— i ■ impossible to derive the corresponding dual theory. Even if a dual field theory can be constructed, it is hard to prove 

' the equivalence of physics at quantum criticality between the original and dual field theories Q ■ 
O'N ■ However, the situation becomes better in one dimension Q. Starting from an interacting one-dimensional model, 
one can perform the weak coupling analysis to find self-consistent equations for self-energy corrections of electrons, 
renormalized interactions, and vertex corrections. Resorting to the Ward identity, one obtains asymptotically exact 
effective interactions, which allow us to find the electron Green's function. It turns out not to have poles but to 
] show two branch cuts, implying the spin-charge separation and identified with Luttinger liquid. This diagrammatic 
analysis can be reformulated in terms of topologically nontrivial excitations called spinons and holons, which carry 
only spin and charge quantum numbers, respectively. This duality transformation is called bosonization, where duality 
equations are highly nonlinear and nonlocal as the case of Ising model. Resorting to such dual variables, one can 
construct a dual field theory in terms of topological excitations, which results in asymptotically the same Green's 
Jv>( , function as that of the diagrammatic approach. 

In this study we address essentially the same issue on non-perturbative dynamics in the Fermi-surface problem. 
Before 2009, it has been believed that the 1/N a expansion, where N a denotes the spin degeneracy to extend from f, 4- 
to 1, 2, N a , is a controllable method in understanding quantum criticality of itinerant electron systems, sometimes 
referred as the Eliashberg framework or a self-consistent RPA (random phase approximation) theory jj]. This belief 
in suspense has been broken in 2009, where all planar diagrams turn out to be in the "same" order even for the 1/ATg- 
expansion This conclusion means that vertex corrections associated with planar diagrams should be introduced 
up to an infinite order in order to obtain reliable critical exponents near quantum criticality Q. 

We investigate two-dimensional antiferromagnetic quantum criticality in itinerant electron systems. Hertz-Moriya- 
Millis theory with dynamical critical exponent z = 2 has been proposed to describe the antiferromagnetic quantum 
criticality. However, we show that the dynamical critical exponent changes from z = 2 to z — 3 at low temperatures, 
which results from effective long-range interactions between spin fluctuations, generated by Fermi-surface fluctuations 
beyond the Eliashberg framework. We argue that the underlying physics for the z = 3 antiferromagnetic quantum 
criticality is the emergence of fermionized skyrmion excitations at low energies, which form a Fermi surface, referred 
as a skyrmion liquid, where the interplay between itinerant electrons and skyrmions is argued to allow fermionized 
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skyrmions. We construct a dual field theory in terms of skyrmion excitations, and obtain the z = 3 antiferromagnetic 
quantum criticality. This demonstration suggests the expected but nontrivial consistency between weak-coupling and 
strong-coupling theories. 

II. HERTZ-MORIYA-MILLIS-CHUBUKOV THEORY REVISITED 
A. Spin-fermion model 

A "standard" model for the spin-density-wave quantum criticality is the spin-fermion model, given by 
Z = J Dc na D<t>e-^ d ^ dd r^ 

C = aL(d T - ii c )c na + -L|(V - iA)c na \ 2 + X^ r cj, ■ c\ a a + 4>(-d* - v 2 V 2 + r 2 )</>, (1) 

where c na is an electron field with spin a =f, \. and n — 1, N, and cf> is a spin-density field with three components. 
fi c is an electron chemical potential, and m c is an electron band mass. A is an electro-magnetic vector potential field, 
externally applied and not considered in this study, is a spin- wave velocity, and £ is a spin-spin correlation length. 
The antiferromagnetic quantum critical point is defined as the renormalized correlation length diverges, i.e., £ — > oo. 
A^ is a coupling constant between spin fluctuations and itinerant electrons, and Q is a nesting wave vector. Although 
we focus on two-dimensional quantum criticality, we leave the dimension to be represented by d. 



B. Eliashberg theory 

The spin-fermion model has been solved in the Eliashberg framework, where only self-energy corrections of electrons 
and spin fluctuations are taken into account self-consistently while vertex corrections are neglected Q ■ Self-consistent 
equations are given by 

E(fe, iu) = 3A|~ G( > k + 1^ + iQ )x(q, 

3A 2 i 

iu k 

where £(fc,iu;) and H(q,ifl) are self-energy corrections of electrons and spin fluctuations, respectively, and G(k,iu>) 
and x(<7,if2) are renormalized electron Green's function and spin susceptibility, given by Dyson equations 

G(k, iuS 



iu: + fj, c — efc — £(fe, iui) ' 

x(9 ' in) = n» + ^ + ^ + n( 9 ,in)' (3) 

respectively. These equations can be solved easily at the antiferromagnetic quantum critical point, solutions of 
which are U(q,iCl) cx NXlvp 1 ]^ for the spin-fluctuation self-energy and E(iw) oc A^A^^^/jaJf for the electron self- 
energy in two dimensions, respectively, where Np is the density of states at the Fermi energy. We note that the 
scaling relation between momentum and frequency in the spin- fluctuation propagator changes from il 2 ~ v 2 ^ |q| 2 to 
|f2| ~ [v 2 v <p j '(A A^,)] |q| 2 at low frequencies due to the spin-fluctuation self-energy. As a result, we obtain the dynamical 
critical exponent of z — 2 for this antiferromagnetic quantum critical point. 



C. Hertz-Moriya-Millis-Chubukov theory 



The z — 2 quantum critical liquid state, described by the Eliashberg theory for spin fluctuations, has been proposed 
to be the zeroth-order mother state at the quantum critical point as the Fermi liquid state away from quantum critical 
points Q. In order to go beyond the Eliashberg framework, we integrate over electronic degrees of freedom, their 
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dynamics of which is given by the Eliashberg theory, and expand the resulting logarithmic term up to the fourth order 
in spin fluctuations. As a result, we obtain an effective field theory for spin fluctuations, 



S *ff = \Y. J ■^<f>r(M,q)(n 2 + vlq 2 +^\n\+C 2 )<f>r(-M,-q) 



1 ^ f d d q I ^ f d d p 1 ^ f d d p' . . 

Wr^J Wy-^J ^M^,q;iv,p^,p) 

% \ L lis %\j 



{[(/>r(iv + itl/2,p+q/2) ■ 4> r (-iv + ifi/2, -p + q/2)][<f> r (iv' - iQ/2,p' - q/2) ■ <j> r (-iv' - iti/2, -p' - q/2)] 

+[4> r (iv + i0,/2,p + q/2) x <f> r {-iv + iCl/2,-p + q/2)} ■ [4> r {iv' - iQ/2,p' - q/2) x <j> r (-iv' - iCl/2, -p' - q/2)}\, 

(4) 

sometimes referred as the Hertz-Moriya-Millis-Chubukov theory. 7^ = NX^v^ 1 is the Landau damping coefficient, 
given by Eq. (2) with Eq. (3), and £ is a renormalized correlation length. Xi(ifl,q;iiy,p;ii>',p') is a renormalized 
interaction vertex with four legs, given by Q 

A4 (ifi, q; iv, p; iv' ', p') 

(i V ( 7^-^G(fc, iu)G{k +p q/2, iu + iv - iQ/2)G(k - q,iw - ifl)G(k - p' - q/2, ico - iv' - iCl/2)\, 

(5) 



f r f 
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where G(k,iuj) is the Eliashberg electron Green's function, given by Eq. (3). 

An essential aspect is that the renormalized vertex of \^(i$l,q;iv,p;iv' ,p') depends on the momentum of q in a 
singular way, giving rise to long range interactions effectively in real space. More precisely, such long range interactions 
occur when the frequency of f2 becomes finite. If one considers the limit of f2 — > 0, this renormalized interaction 
recovers a constant interaction- vertex, referred as the Hertz-Moriya-Millis theory. 

We investigate the role of such nonlocal correlations in non-perturbative dynamics of spin fluctuations, which 
originate from Fermi-surface fluctuations beyond the Eliashberg approximation. Actually, their effects have been 
examined perturbatively, where vertex corrections associated with the nonlocal term are introduced into the Eliashberg 
spin susceptibility 6]. The dynamical critical exponent turns out to get a correction term, resulting in z 7^ 2. 
On the other hand, we perform the RPA analysis beyond the previous study, where the renormalization for the 
A4 (ifl, q; iv, p; iv' ,p') vertex is taken into account. This non-perturbative analysis gives rise to z = 3 antiferromagnetic 
quantum criticality. 

D. Molecular approximation 

In order to renormalize the interaction vertex of X^{iQ.,q;iv,p;iv' ,p'), it is necessary to simplify the nonlocal 
interaction term, but keeping an essential feature. An idea is to introduce "molecular" field variables, which consist 
of composites of original spin variables. Although this procedure can be regarded to be quite conventional, our main 
approximation is to neglect their internal structures and to keep their center-of-mass dynamics only. 

Performing the Fourier transformation for Eq. (4), we start from the Hertz-Moriya-Millis-Chubukov theory in real 
space 



S eff =J*dr J d d rMr,r)[-d 2 T -vlV 2 r +^^dl + r 2 )4>r(r,T) 

I dn / d d ri / dr 2 / d d r 2 / dr 3 / d d r 3 / dr 4 / d d r 4 
J Jo J Jo J Jo J 



A 4 [?"2 - r 1 ,r 2 - Ti ; r 4 - r 3 , r 4 - r 3 ; (r 4 + r 3 )/2 - (r 2 + r ± )/2, (74 + t 3 )/2 - (r 2 + t x )/2] 

|[</»r(n,n) • 4>r{r2,T 2 )][ct>r{r- i ,T 3 ,) ■ </> r (r 4 ,T 4 )] + [0 r (n,Ti) X r (r 2 ,T 2 )] ■ [(f> r (r 3 ,T 3 ) x r (r 4 ,r 4 )]|. (6) 



It is straightforward to rewrite this field theory in terms of center-of-mass and relative coordinates, given by 

S ef f = f dr J d d r<j> r {r,T)(-d 2 T -vlV 2 r +^^& T +C 2 )<t>r(r,T) 

rP r rP r rP r rP 

dn / d d n / dr 2 / d d r 2 / dTx \ d d R 1 / 
>o J Jo J Jo J Jo 



+ [ dn [ d*n [ dr 2 f d d r 2 f dTifd d R 1 f dT 2 I ' d d R 2 \ A [r l ,T l -r 2 ,T 2 ;R 2 - R U T 2 -Tx\ 
Jo J Jo J Jo J Jo J 

{[0 r (fli +ri/2,Ti +n/2) • MRi - n/2, 71 - ri/2)] [ct> r {R 2 + r 2 /2,T 2 + r 2 /2) • 4> r (R 2 - r 2 /2, T 2 - r 2 /2)] 

+ [<j> r (Ri + n/2, Ti + ti/2) x cj )r (R 1 - n/2, 71 - n/2)] • [0 r (iJ 2 + r 2 /2, T 2 + r 2 /2) x r (J? 2 - r 2 /2, T 2 - r 2 /2)^) 

where R\ = ri + V2 , 7^ = Tl + T2 and H 2 = r3 + Ti ; 7^ = T3 + T4 are center-of-mass coordinates and r 2 — n — ► n, 
t 2 — n — >■ ti and r 4 — r 3 — ► r 2 , r 4 — t 3 — ► r 2 are relative coordinates. 
We introduce molecular field variables, 



<t> r {Ri + n/2, 71 + n/2) • <M-Ri - n/2, 71 - n/2) \ — > V(n, n; Hi, 71), 

r (Hi+n/2,Ti+ri/2) x -n/2,Ti-ri/2)\ — > Ar(n,n;-Ri,Ti), (8) 



where V(ri, n; 7i) represents valence bond fluctuations, and X(n,Tx;R\,Ti) expresses vector spin-chirality fluc- 
tuations. 

Performing the Hubbard-Stratonovich transformation and doing the Fourier transformation, we obtain 

S *ff = E E ^ + V W + + ^~ q ' 

if! q 

+EE^EEt;EE^EE M-«n> -(« + <*')A + ^oa 

iQ. q P iQ' q' P iOi qi P jfi'j qj 

{^V(qi,iOi; <^,ini)V((<| + <?')A (ifi + ifi')/2; -qi, -iSl[) - iV(«n,ifii;g - g'.ifi - ift')[0r(<?', ■ 0r(-<7, -ifi)] 

+ jAT(gi,ifii; gi.ini) • AT((cj> + q')/2, (in + HI')/ 2; -q[, -in[) - iAT^.ifii; q - q' , - ifi') • [0 r (q', ift') x r (-q, -iQ)]j. 

(9) 

Up to now, all procedures from Eq. (4) to Eq. (9) are exact. 

An essential approximation is to neglect internal dynamics of molecular field variables and to keep their center-of- 
mass dynamics. Actually, this approximation has been rather conventionally utilized because internal dynamics can 
be regarded as fast degrees of freedom. However, we confess that the validity of this approximation is not proven, 
particulary, at the quantum critical point. If sizes of molecules become enhanced and the probability for overlapping 
between them gets larger at the quantum critical point, our approximation to view such molecules as point particles 
will not be valid any more. 

Performing the point-particle approximation for such composite field variables, we obtain 



S *ff ~ E E Mq, in ) (° 2 + V W + 7V>M + £ 2 ) M-q, 



in q 



+ EE{^%^ V (<?, in)v(- q , -m) + ^Rx( g , m) ■ x(- g , -m) 

iU q 

-*(tj E E A 4(9i' ^i)) !EE V («- ifi - in'Mrtf, m') ■ <j> r {- q , -in)] 

ifi'j q[ iQ' q' 

- l G E E M^i)) 4 E E - in - in') ■ [<j> r { q ', in') x <m- 9 , (10) 

P in[ q[ P in' q' 

Rescaling V(q - q',in - in') — -> V(g - <?',ift - *fi')/(p£ini A 4 (q 1 ,iO / 1 )), - - iQ') — > ^(q - 

g'.in-ifi')/^ £\ n , A 4 (qi, ifii)), and A 4 (q, ifi) — >• ^ J2tn[ J2 q [ ^(Qi^^ii) A 4 (q, ifi), we reach the following 
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expression 



z eff = J Dcj) r (q,in)DV{q,in)DX{q,in)exp{- ps eff [4> r (q : in),v(q,in),x{q,in)}y 
s -ff = E E (n 2 + vj q 2 + ^\n\ + r 2 ) 0r(-g, -in) 



iO q 
if! q 

-^^^V(q-q',iO-ffl')[0r(g',^') • 0r(-<Z, -in)] - i„ E E ^ ~ g '' m ~ ' x 0r(-g, -in)]}. 

iH' q' zQ' q' 

(ii) 

One can evaluate the four-point interaction vertex, given by 0] 

mi 

A 4 ( 9 , l O)« A 4 ^ (12) 

asymptotically, which implies that this interaction is marginal in two dimensions. 

Integrating over molecular fields and performing the Fourier transformation in Eq. (11), it is not difficult to 
understand the reason why we are saying that only the center-of-mass dynamics is considered, 

S eff =J dr J ^r0 r (r,T)(-a2-^ V 2 +7V) yra2 +r 2 )0 r (r,r) 

+ / dn [ d d ri [ dr 2 I ' d d r 2 [ dT 3 [d d r 3 \— r[0 r (ri,n) • <j>r{r^n)][4> r {r 2 ,T2) ■ 4>r(ri,n)] 

Jo J Jo J Jo J L A4(ri - r 2 ,ri - r 2 ) 

+ T-? A(j>r(ri,n) x r (r 3 ,r 3 )] • [0 r (r 2 ,r 2 ) x r (r 3 ,r 3 )]j. (13) 

M{ri - r 2 ,n - r 2 ) J 

An interesting point is that the interaction vertex appears in the denominator instead of the numerator. Of course, 
this is consistent with its dimension, considering how it rescales. This relation reminds us of the weak-coupling and 
strong-coupling duality. 

III. BEYOND THE ELIASHBERG FRAMEWORK 
A. Renormalization for the interaction vertex in the RPA level 

An important point beyond all previous studies is to introduce renormalization of the interaction vertex. We 
perform the renormalization within the RPA analysis. It is straightforward to see that no renormalization occurs 
from vector spin-chirality fluctuations, given by 

iVL q iQ,' q' iQ," q" iQ 1 " q'" 

X(q - q', in - iSl')X(q" - q'",iQ" - in'") = 0, (14) 

where the subscript "c" denotes "connected" for the diagrammatic expansion. On the other hand, valence bond 
fluctuations give rise to 

s lt « \ E E \ E E \ E E \ E E(^(«'> i0 ') ■ m-q, -m [M*"', • m-q", -m')] 

zSl q iQ' q' iQ" q" iQ"' q'" 

V(q - q', in - in')v(q" - q'", in" - in'") 

~ E E g E E(^( g ' ifi ) ' (Mi 1 , • M-q', v (<? - ifi - »fi')v(-g + <?', -in + in') 

iQ q iQ' q' 

= \ E E{ 7? E E w(9 + in + ifi ')^(<z', »n')}v(«, in)V(-g, -*n), (15) 

in g iQ' q' 

where ^(g, in) = ^f2 2 + v^q 2 + 7^|n|^ is the spin-fluctuation propagator. 



6 



B. Self-energy corrections to antiferromagnetic fluctuations 



The renormalized interaction vertex causes the self-energy correction for spin fluctuations, given by 

S L ~ 2 E E \ EE( V (« - <*', in - ifi')v(-q + q', -in + in')) ^(<t> r (q' , in') • <M-<z', -in') 



i£7 q 



iH' q' 



4> r {q, in) • r (- g , -<n) = E E{^ E E G ^(9 - *n - in')^*?', in')}<M<z, in) ■ <M-<z, -in), (16) 



if) q i£"2' q' 

where the RPA propagator for valence bond fluctuations is 
G v (q,iQ) 



1 



+ H ^(9 + *n + in")M<z", 

As a result, the effective action for spin fluctuations reads 



(17) 



S tff =^2^2<l>r(q,ify 

zfi q 

1 1 



n 2 + vW 



■7vin| 



+ 



2/? 



y y 

/ j Z ^ A4(q — 



g^{q',iQ!) 



q' 



q' .iil — iil' ) 



4 d 



£„» 5^(9 + q", ^ - in' + iQ")9<t>(<l" , 



r (-g,-in). (18) 



This expression shows the main point of our study. If the RPA renormalization is not introduced, the self-energy 
correction for spin fluctuations is reduced to that of the previous study @. Indeed, one can find logarithmic sin- 
gularities from the "bare" interaction. However, we find that the self-energy correction in the RPA valence-bond 
fluctuation propagator overcomes the "bare" part. In other words, we obtain Gy(q,iQ) ps {j-^J2m" J2 q " d^iQ + 

q",«n + ifl'^g^q" , in")} . Inserting this expression into Eq. (18), we find 



c4> 



EE^(<?,*n) 7v ,|n| + ^g 2 - 

ifl q 



H^ 2 

c rig* 



ln(— iO) <fi r (—q,—i£l). 



(19) 



Detailed calculations are shown in appendix. As shown in this expression, z — 3 antiferromagnetic quantum criticality 
emerges in q < q c , where q c is a characteristic momentum given by the comparison between the first term and others. 
This leads us to conclude that the z = 3 Hertz-Moriya-Millis-Chubukov theory describes the antiferromagnetic 
quantum critical point at low temperatures. 



IV. A STRONG-COUPLING APPROACH 



A. An effective field theory 



It is not easy to understand the underlying physics for the emergence of the z = 3 antiferromagnetic quantum 
criticality from the z = 2 Hertz-Moriya-Millis-Chubukov theory. Here, "to understand the underlying physics" means 
to identify elementary excitations, resulting in z — 3. As described before, the renormalization of the interaction 
vertex from valence bond fluctuations is the source of the z = 3 quantum criticality. In this respect our question is 
as follows. Can we find the dual field variable of valence bond fluctuations? We claim that fermionized skyrmion 
excitations play the role of valence bond fluctuations. We show that fermionized skyrmion excitations with their 
Fermi surface result in z — 3 quantum criticality. 

We start from the following spin-fermion model without the kinetic energy term of spin fluctuations, 

Z=[Dc n<7 Dcf ) e-fo^rf^rC > 

C = c) na (d T - n c )c nrJ + 7^-|(V - iA)c niJ \ 2 + \ 4> e lQ r 4> ■ c\ a a a pc n0 , (20) 

regarded as the strong coupling limit of Eq. (1). We introduce the CP 1 representation for the angular field variable, 
where the amplitude field is assumed to be frozen in the low energy limit, given by [8J 

\^ Q r cf> ■ a afj = \^ Q - r \<j>\U ai a z lS Ul„. (21) 
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U = ( ^ ^ ] is an SU(2) matrix field to represent the angular fluctuation, where — e 1 * cos § and zj_ = 

e % ~5 sin | denote complex boson fields, identified with spinons. 

An essential ansatz in the strong coupling limit of is that the spin of an electron field is screened from the 
angular fluctuation of the spin-fluctuation field, represented as follows 



i/w = Ul ff c n p. (22) 



ip na is a holon field which carries the same charge quantum number as an electron while it does not carry the spin 
quantum number. Based on this ansatz, we rewrite the spin-fermion model in the strong coupling limit as follows 



+X 4> e iQ - r \cj>\cx^ na + ^-^(d r U^ + [tfdrUpA (d T U + U[(d r U^)U}) </W - ^-d T (^M, (23) 
2m c \ /07 V / 7/3 zm c 

where holons and spinons appear to realize the spin-charge separation. Here, we omit the electromagnetic vector 
potential for the time being. We emphasize that this expression is just a reformulation of Eq. (20), regarded as the 
change of variables. 

Introducing the nonabelian gauge field of A v a g = —i\{d v U^}U] a p, one can rewrite the above as follows 
Z = J D$ w DU a pDJ^ p 8{Ul 1 Utf-8 a e)6(J^p + i[^^^ 

+ 2^L (drU^ - iA r aS ul) (d r U lP + iU 7 sAlp)i> n0 , (24) 

where SU(2) gauge symmetry becomes shown explicitly. 

Although exactness is still guaranteed in Eq. (24), it is rather complicated to handle this effective field theory. In 
this respect we perform the U(l) projection, considering the subset of the SU(2) gauge theory, given by 

Z = J Dij na Dz a Da^(\z a \ 2 -l)S(d r a r )e-^ d -I d2 ^ 7 C = £^ + £ z , 

= iptvidr - Mr + A0e lQ " r |^|cr - iaa T )ip na + -— \(d r - i<ra r - iA r )ip nrT \ 2 , 

2m c 

= T-I(0r - ™ T )z a \ 2 + 1 EL.\(a r - ia r )z a \ 2 (25) 
2u c Zm c 

with \x r = \i c + m^ 1 . One may understand this procedure as the staggered-flux ansatz in the SU(2) slave-boson 
theory where the SU(2) gauge symmetry is reduced to the U(l) symmetry. It is interesting to see that the internal 
gauge field a r couples to the spin current of the renormalized electron field while the electromagnetic field A r does to 
the charge current. The spinon part is reduced to the CP 1 representation of the 0(3) nonlinear a model @ with a 
stiffness parameter of p s = (cn[)J la .il; ncr ) , where the time derivative term is added explicitly. This time derivative term 
is expected to appear from quantum corrections, i.e., the self-energy correction to the spinon dynamics. 



B. Duality transformation 

To pull out dynamics of skyrmions, we perform the duality transformation. Unfortunately, the duality transforma- 
tion for the SU(2) case is not known. We take the easy plane approximation 

z a = ±e^. (26) 

In this representation the topological excitation corresponding to a vortex is a meron instead of a skyrmion [l(| ■ Since 
a meron solution can be regarded as a half skyrmion, we will introduce a meron pair as a skyrmion [llj . 
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Resorting to the easy plane approximation, we obtain 

£ = i^nA^r - Mr + \^ lQ ' r \<f>\<j - iaa T )ip nlT + -^—\(d r - iaa r - iA T )ip na \ 2 + -^—(d r (j> a - a T ) 2 + -^-\d r <j> a - a r \ 2 . 

(27) 

Performing the duality transformation, we obtain the dual Lagrangian 

C D = ^ nrT {d T - fi r + A0e jQ ' r |0|cr - i<ra T )ip na + ^~l(^r - iaa r - iA r )ip na \ 2 

+ y (e^^O + — ^(e rAlI ,<9 M c£) 2 - iJ a Dll cl - ia^xducl, (28) 

where Jp is a meron current and is the dual gauge field representing spin wave excitations. 
Introducing the second quantization for meron dynamics (lpj . we obtain 

C D = ipl a {d T - \± r + X^e %Q ' r \<f}\i7 - iaa T )ip na + ^HO^r - iva T - iA r )ip na \ 2 

+ \(d, ~ ic°W D \ 2 + m 2 D \<S>° D \ 2 + fW D \ 2 ) 2 + V an [*° D ] + y (e.^^O 2 + ^(e r ^d^f, (29) 

where $J> is a meron field and \4 n [$£)] is an anisotropy potential for merons. 

To find skyrmion dynamics, we rewrite the above meron Lagrangian in terms of meron pairs. This derivation is far 
from rigorous, but sometimes utilized in order to make physical excitations 11 [. Considering — > e l0<T , the meron 
kinetic energy term is 

[d,e a - C p 2 = (d,9 t - cl) 2 + {8,6, - ci) 2 

= + W H + ci]) 2 + l([d,e t - 8,6,] - [cl - c^]) 2 . (30) 

The spin wave term is rewritten as follows 

2 \^TfiuOfj,C u ) + \.^rnvOft,C v ) 



.( e dc^ + e d c l Y + —(e d c T - e d cA' 
'-(e d c T + e d c l Y + — 

4p s \ / 4p s 



Uc 
4 

IUr .( f dA + f d c l Y + —(e dA-e 8 r l Y CUl 



As a result, we reach 



where 



+ 1(5, - ic+)*+| a + m+ 2 |<I>+| 2 + ^|d>+| 4 + |( £t ^c+) 2 + g^d^) 2 

+ l(^ - ^c-)<5> D \ 2 + m D 2 \<S> D \ 2 + ^ D \ 4 + ^(e T ^9 M c") 2 + g^^^c") 2 + V an [^>+, (32) 

*± = — ^ 4 = JL ^f L - (33) 

It is interesting to see that the field docs not couple to dynamics of renormalized electrons ij) a directly. Only the 
field couples to fermions via the mutual Chern-Simons term. Integrating over c^, and a M fields, we will obtain 
an effective field theory for renormalized electrons interacting with skyrmions. It is natural to expect the emergence 
of nonlocal interactions between spin currents of renormalized electrons (holons) and skyrmion currents. Integrating 
over holon excitations, we reach a dual filed theory for skyrmion dynamics. Unfortunately, this procedure cannot be 
performed rigorously. We propose the following skyrmion field theory 

S sk = J drj d 2 r{ Msfe $t(d T _ JCt )$ s + _ iCM )$ s |2 + m 2|$ s |2 + ^| $a |4 + ^^g^f + ^(^a^) 2 }, 

(34) 
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where $^ i s replaced with $ s and the — superscript in c~ is omitted. 

It is our proposal to introduce the linear-time derivative term into Eq. (34). Then, the linear-time derivative 
term becomes more relevant than the second-time derivative term at low energies. Unfortunately, we cannot prove 
the emergence of the Galileian invariance (the linear-time derivative) from the relativistic invariance (the second-time 
derivative) at present. In this respect one can call the emergence of the Galileian invariance in the presence of itinerant 
electrons as our speculation. 

There must be an underlying physical mechanism for this linear-time derivative term, which breaks the particle-hole 
symmetry for skyrmion excitations. It has been demonstrated that the density of skyrmion excitations is finite at an 
antiferromagnetic quantum critical point without itinerant electrons, i.e. in an insulating antiferromagnet [l2j . Of 
course, only the second-time derivative term is allowed in this case, i.e., particle- hole symmetric, which means that an 
equal number of skyrmion and anti-skyrmion excitations exists at this quantum critical point. Our problem is what 
happens on the particle-hole symmetry in skyrmion excitations at the quantum critical point where itinerant charge 
carriers are introduced. This is a long-standing problem, where non-perturbative effects from interactions between 
itinerant electrons and "many" topologically nontrivial excitations should be taken into account on equal footing. 
Frankly speaking, we do not have any reliable mathematical tools for the description of such interactions. 

Our speculation is that the presence of itinerant electrons will induce the particle-hole symmetry breaking in 
the skyrmion sector because itinerant electrons favor skyrmion excitations (or anti-skyrmions, i.e., one of the two). 
The physical mechanism is as follows. When itinerant electrons move in the background of skyrmions, they feel an 
effective magnetic flux, which quenches the kinetic energy of electrons. Our expectation is that the gain in the electron 
kinetic energy contribution can overcome the energy cost due to the particlc-holc symmetry breaking in the skyrmion 
sector. As a result, the particle-hole symmetry breaking is favorable in the respect of the total energy. Actually, this 
mechanism has been realized in the system of frustrated magnets, where the presence of itinerant electrons leads the 
co-planar ordering in the triangular antiferromagnet to be ordered into an out-of-plane way, which corresponds to a 
spin chiral order [13j . 

Possibility of statistical transmutation for vortices has been discussed in the vortex liquid phase [l4j and at quantum 
criticality in geometrically frustrated spin systems [TBI - Fermionization of skyrmion excitations can be performed in 
the same way as that of vortices. The key point is that the Chern-Simons term becomes irrelevant at quantum 
criticality [lj, , implying that the quantum statistics for skyrmions or vortices may not be well defined at quantum 
criticality. Actually, such excitations are strongly interacting at criticality, where we are not allowed to pin down 
elementary excitations clearly. We reach an effective field theory of fermionic skyrmions, 

r fir r /- j ^ tti "l 

S eff = dr d 2 r\ ipl(d T - Ity- ic T )ip s + |(V - ic)ijj s \ 2 + -£-{e T u, v du,c v ) 2 + TT^i^^d^) 2 \, (35) 

Jo J 1 zm sk Z Zp s > 

where ip s represents the fermionic skyrmion field with the skyrmion chemical potential jj,^ and the band mass m s k oc 
Hsk- We would like to emphasize that Eq. (35) can be derived from Eq. (34) at least formally via the Chern-Simons 
transformation [l4| [l5| . 

Another important point with the particle-hole symmetry breaking term is that the skyrmion chemical potential 
is assumed to be finite. As discussed before, the skyrmion density has been shown to be finite in the 0(3) nonlinear 
a model without itinerant electrons 12]. It is natural to expect that this will hold even in the presence of itinerant 
electrons. Then, Landau damping for gauge fluctuations emerges from particle-hole excitations of fermionic skyrmions 
near their Fermi surface, resulting in the z — 3 dynamics for spin fluctuations. 

Although the z = 3 antiferromagnetic quantum criticality appears in the fermionized skyrmion ansatz with the 
skyrmion Fermi surface, it is still not clear how this z = 3 quantum criticality can be connected with that of the 
previous diagrammatic approach. It has been shown that the skyrmion field has the same symmetry as the valence 
bond field in an insulating antiferromagnet 

[HUE 

[l7j . Resorting to the fermion representation for the spin operator, 
one can construct an effective Dirac theory with an enhanced emergent symmetry, compared with the Heisenberg 
model with the 0(3) spin symmetry [l6l - [T8j . Based on the Dirac theory, one can construct the valence bond operator 
and skyrmion operator in terms of the Dirac spinor. Since the symmetry property is well defined for Dirac fermions, 
it is straightforward to investigate symmetries of both valence bond and skyrmion operators. They turn out to be 
the same as each other. One can perform the same work in the boson representation for spin, where the role of the 
spin Berry phase is crucial to assign the valence bond quantum number to a skyrmion [10(. On the other hand, it 
is not verified yet whether skyrmion excitations can be identified with valence bond fluctuations in the presence of 
itinerant electrons. In particular, the role of Berry phase is not clarified yet. However, it is expected that valence bond 
fluctuations will be deeply related with skyrmion fluctuations. Actually, the skyrmion operator has been constructed 
in terms of fermion bilinear operators, identified with the valence bond operator in itinerant antiferromagnets [l9| . 
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V. DISCUSSION AND SUMMARY 

In summary, we showed that the dynamical critical exponent can change from z — 2 to z = 3 at low temperatures in 
antiferromagnetic quantum criticality, where nonlocal interactions between spin fluctuations, which result from Fermi 
surface fluctuations, are responsible for the change of the dynamics of spin-fluctuations. In particular, we claimed 
that renormalization in the nonlocal interaction vertex should be introduced because the renormalized part turns 
out to be larger than the bare (unrenormalized) interaction. We performed the renormalization in the RPA level, 
where valence bond fluctuations are introduced explicitly and the nonlocal interaction vertex is identified with a bare 
propagator of valence bond fluctuations. As a result, the RPA renormalization in valence bond fluctuations gives rise 
to the z = 3 antiferromagnetic quantum criticality. 

We tried to understand the underlying physics for the emergence of z — 3 antiferromagnetic quantum criticality. 
Hinted from the fact that valence bond fluctuations can be identified with skyrmion excitations, we constructed a 
dual field theory in terms of skyrmion excitations. The key approximation is to introduce the particle-hole symmetry 
breaking term, i.e., the linear-time derivative term instead of the second-order time-derivative term. This construction 
is based on the fact that the interplay between itinerant electrons and skyrmion excitations may favor skyrmion exci- 
tations more instead of anti-skyrmion fluctuations. Allowing the statistical transmutation for skyrmion fluctuations, 
where the Chern-Simons term turns out to be irrelevant at the quantum critical point, we found a dual theory in 
terms of fermionized skyrmion excitations interacting with spin- wave fluctuations. If the density of skyrmions is finite 
at the quantum critical point, fermionized skyrmions will form their Fermi surface. As a result, Landau damping 
dynamics for spin-wave excitations is introduced to cause the z = 3 antiferromagnetic quantum criticality. 

We propose that this z = 2 to z — 3 crossover in two dimensions may explain the crossover behavior in the specific 
heat measurement of YbRIi2Si2, where the logarithmic divergent behavior of the specific heat coefficient turns into 
the power-law divergent behavior of T~ 2 / 3 at low temperatures [2p|. |2lj. 

An important missing point in our study is to neglect higher-order terms beyond the fourth order expansion. Such 
higher-order terms have been argued to be important because the two-dimensional z — 2 spin-fluctuation theory 
turns out to allow infinite number of marginal interactions, which originate from singular momentum and frequency 
dependencies for interaction vertices. Based on this observation, it was claimed that an anomalous dimension for spin 
fluctuations appears beyond the Eliashberg framework (22[. We do not exclude this possibility. Actually, the spin 

susceptibility can be Xr(<Z>&n) = (^7i>l^l + v2 ,Q 2 — ln(— instead of Xr(Q,i^) — (t^I^I + v2 Q 2 

—f- 4 4 ln(— iW in Eq. (19), where 7 is the anomalous dimension of spin fluctuations. We leave the study beyond 
the fourth order as a future work. 
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Appendix: Self-energy corrections to antiferromagnetic fluctuations 

The RPA self-energy correction for valence bond fluctuations is given by 



\ E E 94>(q -1' + 1", in - & + in") g ^(q", in" 



ft iU" q" 



ivy I I 

ft W ^ (n - n' + tt") 2 + vf(q -q' + q") 2 + 7^ - Q> + Q"\ ^" 2 + ^l" 2 + 7^"\ 

^ ^ poo p2ix ^ ^ 

Wfl^X dq " q "Jo d9 vlq" 2 + 2vy\q - q'\ cos6 + v^q - q'f + 7^ - fi' + Q"\ v%q"* + 

In ^ (AT) 



47t 7 ^ n - n> 

where we resorted to q sa q' in the last line. C is a positive numerical constant. This renormalized contribution turns 
out to be much larger than the bare vertex of \i(q, ifl) because Jj^lr -C 1 should be satisfied in the expansion. As a 
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result, the self-energy correction in spin fluctuations is given by 



- w 



g^{q',iQ!) 



Pj^^f-p E,» - g' + <?", in - i«' + <fi")$*(g", *«") 

94>{q,i^') w 4?r 7^ 1 \ ^ 1_ 



^ ?E;n» E g » <^(<?",*n - in' + in'Os^ifi") c /3^ ln (_^_^ ^g 2 + 7v ,|fi'| 

= gf /l nH 0)ln( -^ 2 4 + ^ g4 ) - p dy W 2t -^^) -P^(2,-^» )}. (A.2) 

Considering that the PolyLog function is not singular at low momenta and frequencies, we expand the first log term 
in |fi|/|q| 2 <C 1 and find the expression of Eq. (19). 
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